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In the search of new electrocaloric materials: Fast ion conductors 
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We analyse the effects of applying an electric field on the critical temperature, Ts, at which superi- 
onicity appears in archetypal fast ion conductor CaF 2 by means of molecular dynamics simulations. 
We find that the onset of superionicity can be reduced by about 100 K when relatively small electric 
fields of ~ 50 KV cm“^ are employed. Under large enough electric fields, however, ionic conductivity 
is depleted. The normal to superionic phase transition is characterised by a large increase of entropy, 
thereby sizeable electrocaloric effects can be realised in fast ion conductors that are promising for 
solid-state cooling applications. 


Introduction. In electrocaloric materials, the adiabatic 
switch of an electric field causes a change in the temper¬ 
ature of the system that is equal to 



where Cp represents the heat capacity, S the entropy, 
and E the external electric field. Ferroelectrics are the 
archetypal electrocaloric compounds>i They exhibit an 
spontaneous electrical polarization below a certain crit¬ 
ical temperature, Tc, that can be shifted by an exter¬ 
nal electric field. Electrocaloric effects are observed at 
T > Tc, when the system is paraelectric but responds to 
the presence of an external electric field by aligning its 
dipole moments parallel. In those conditions, the adia¬ 
batic switch of an electric field causes a positive AT in 
the crystal because the entropy in the ordered state is 
smaller than in the paraelectric phase [AS' < 0 hence 
AT > 0, see Eq. (HI)]- Conversely, when the electric field 
is adiabatically removed the material gets cooler (AS > 0 
hence AT < 0). Several thermal cycles based on the elec¬ 
trocaloric effect just described can be engineered, which 
are promising for solid-state refrigeration applications.— 
Nevertheless, ferroelectric materials present a series of 
technical issues, like for instance the presence of ferro¬ 
electric domains and leakage currents, that are hindering 
the development of practical electrocaloric applications— 
Therefore, is very desirable to find new electrocaloric ma¬ 
terials with improved cyclability and electrical resistiv¬ 
ity features. Promising materials rivaling ferroelectrics, 
overall, must display a large change of entropy upon the 
application of an external electric field. 

Long time ago, Kharkats et al. proposed, based en¬ 
tirely on theoretical arguments, that the application of 
an external electric field could drastically lower the crit¬ 
ical temperature, Tg, at which superionicity appears in 
fast ion conductors (FIC)— Superionicity refers to the 
unusually large mobility (of the order of I 17“^ cm~^) 
that a particular atomic species in a multicomponent 
ionic material acquires when temperature is raised (be¬ 
low its melting point). The normal to superionic phase 
transition is experimentally characterised by a sudden 
increase in the heat capacity and entropy of the crystal 
(A^ ^ 10 J mol“^K“^)— The possibility of externally es- 



FIG. 1. Schematic diagram of an electrocaloric cooling cycle 
involving FIG at T < Ts- “N” and “S” represent the normal 
and superionic states, (a) (b) An electric field is adiabati¬ 

cally applied that triggers superionicity, hence the crystal gets 
cooler, (b) ^ (c) The system receives heat from the environ¬ 
ment, hence T and S increase, (c) —>■ (d) The electric field is 
adiabatically removed, hence the system becomes normal and 
T increases, (d) ^ (a) Heat is ejected from the system. 


timulating superionicity by means of an electric field has 
been, to the best of our knowledge, totally overlooked 
to date. This idea, however, has the potential to mo¬ 
tivate original searches on new electrocaloric materials 
for solid-state cooling applications, as we schematically 
show in Fig. |TJ In this molecular dynamics work, we 
show that Tg can be modified appreciably in archetypal 
fast ion conductor CaF 2 by applying an external electric 
bias. In particular, we show that when E is relatively 
small (^ 10 KV cm“^) Tg can be shifted down by about 
100 K. However, under the action of larger electric fields 
ionic conductivity is depleted. 

Methods. Our {N,P,T) simulations were performed 
with the LAMMPS code, 0 keeping the temperature 
and pressure fluctuating around a set-point value by us¬ 
ing Nose-Hoover thermostats and barostats. Large boxes 
containing 6,144 atoms were simulated over long times of 
~ 100 ps, and periodic boundary conditions were applied 
along the three Cartesian directions. Newton’s equations 
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FIG. 2. Anionic MSD function calculated in CaF 2 considering 
different temperatures and values of the electric field. 


of motion were integrated using the customary Verlet’s 
algorithm and a time-step length of 10“^ ps. A particle- 
particle particle-mesh fc-space solver was used to com¬ 
pute long-range van der Waals and Coulomb interac¬ 
tions and forces beyond a cut-off distance of 12 A at 
each time step. The interactions between ions were mod¬ 
eled with the Born-Mayer-Huggins ion-rigid potential de¬ 
scribed in work Q. The suitability of this approach for 
studying the energy, structural, and superionic proper¬ 
ties of CaF 2 at normal pressure has already been demon- 
stratedi^i^ We simulated the effect of applying an uni¬ 
form external electric field by adding a force equal to 
—gE on each ion (where q is the corresponding charge). 
Following previous works)^^— we identified the onset of 
superionicity with the appearance of a non-zero slope in 
the mean squared displacement function (MSD) calcu¬ 
lated for the fluorine ions. The MSD function is defined 
as {Ar^(t)) = ({ri{t + to) — ri{to))'^) , where ri{t) is the 
position of atom i at time t, to an arbitrary time origin, 
and (• • •) denotes average over F” ions and time origins. 

Results. Figure [2] shows the fluorine MSD function cal¬ 
culated at T < Tg and T > Tg [where Tg = 1350 (25) K 
is the transition temperature computed at E = 0]^ con¬ 
sidering different values of the electric field. E was ap¬ 
plied along the three inequivalent crystalline directions 
[100], [110], [111], and in all the cases we obtained equiv¬ 
alent results. In our analysis we only consider electric 


fields smaller than 500 KV cm~^ since otherwise we found 
that the Ca^+ cations started to drift. At T = 1250 K, 
that is, about 100 K below Tg, the onset of superion¬ 
icity appears at if = 50 KV cm“^. When the module 
of the electric field is further increased, however, ionic 
conductivity disappears. Such a counterintuitive effect 
can be understood as follows. At small values of E, the 
premature creation of Frenkel defects (i.e., simultaneous 
formation of an interstitial-vacancy pair) is estimulated 
by the presence of the external electric field and the ac¬ 
companying increase in the entropy, which minimises the 
free energy of the sytem. Importantly, at these condi¬ 
tions the diffusion of the F“ ions is Brownian (that is, 
(r(t)) = 0), exactly as it is observed in the crystal at 
T > Tg in the absence of an electric bias. Eventually, 
as E is increased, the diffusion of the anions starts be¬ 
ing non-erratic {{r{t)) A 0) and thus the entropy of the 
system decreases and the Coulomb interactions between 
ions are not optimal. Consequently, the anions prefer to 
remain in the ordered normal state, in which the total 
electrostatic energy is most favourable, in order to lower 
the free energy of the crystal. At T = 1450 K, superionic¬ 
ity is fully developed and the effect of applying an electric 
field is to further promote the diffusion of anions. The 
ionic conductivity is maximum at E = 100 KV cm“^, 
however under increasing electric bias the slope of the 
MSD function is reduced (although this is always positive 
and larger than the obtained at E = 0). The origins of 
this effect can be understood in terms of similar entropy 
and electrostatic energy arguments than explained above. 
Finally, we calculated the heat capacity and change of en¬ 
tropy associated to the superionic transition in CaF 2 at 
E < 100 KV cm-i and T < 1350 K. We found that 
the resulting change of temperature assuming adiabatic 
conditions generally is |Ar/r| « 1 %. 

Summary. In conclusion, our molecular dynamics work 
shows that is possible to vary Tg in FIG by applying 
relatively small external electrical fields. This effect has 
the potential to change our paradigm in the search of new 
electrocaloric materials, which are promising for solid- 
state cooling applications. Analogous studies performed 
on similar FIG with lower Tg’s (e.g., Q;-PbF 2 ) are highly 
desirable. 
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